library(tidyverse)
library(sf)
library(terra)
library(openxlsx)
library(mapview)
library(RColorBrewer)
library(osrm)
library(movecost)
library(fs)
library(stringr)
library(synthdid)
mapview::mapviewOptions(fgb=FALSE)
load("data/china.RData")

衛星夜間光データとは

日本の1992年と2018年における夜間光

ntl1992<-terra::rast(x="data/9828827/Harmonized_DN_NTL_1992_calDMSP.tif")
ntl2018<-terra::rast(x="data/9828827/Harmonized_DN_NTL_2018_simVIIRS.tif")
ntl1992 <- ntl1992 %>%
  terra::crop(y=ext(c(123.1,146.9,24.0,45.7)))
ntl2018 <- ntl2018 %>%
  terra::crop(y=ext(c(123.1,146.9,24.0,45.7)))

1992年

terra::plot(ntl1992)

2018年

terra::plot(ntl2018)

Synthetic DiD法とは

\[ \sum \]

北海道拓殖銀行破綻と道内経済への効果検証

jpn1<-sf::read_sf(dsn="data/JPN_shp/gadm41_JPN_2.shp")
mapview::mapview(x=jpn1,
                 zcol="NAME_2",
                 color="transparent")

中国青海-西蔵鉄道開通による経済効果分析